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abstract 

Multiresolution approximation (MRA) of the vector fields on is studied. We intro- 
duced in the Fourier space a triad of vector fields called helical vectors which derived 
from the spherical coordinate system basis. Utilizing the helical vectors, we proved the 
orthogonal decomposition of L^(T^) which is a synthesis of the Hodge decomposition 
of the differential 1- or 2-form on T'^ and the Beltrami decomposition that decompose 
the space of solenoidal vector fields into the eigenspaces of curl operator. In the course 
of proof, a general construction procedure of the divergence-free orthonormal complete 
basis from the basis of scalar function space is presented. Applying this procedure to 
MRA of L^(T'^), we discussed the MRA of vector fields on and the analyticity and 
regularity of vector wavelets. It is conjectured that the solenoidal wavelet basis must 
break r-regular condition, i.e. some wavelet functions cannot be rapidly decreasing func- 
tion because of the inevitable singularities of helical vectors. The localization property 
and spatial structure of solenoidal wavelets derived from the Littlewood-Paley type MRA 
(Meyer's wavelet) are also investigated numerically. 



1 



1 Introduction 



Divergence free vector fields with coherent structures are ubiquitous in a lot of natural 
phenomena, for example, coronal flare of the Sun, dipolar magnetic field of the Earth, or 
the coherent vortices in sufficiently subsonic motions of fluid, for example Great Red Spot 
in Jupiter. Wavelet analysis has been regarded one of the promising tools for surveying 
such coherent structures. Because they are remarkably localized not only in physical 
space but also in Fourier space within the limit of the uncertainty principle. Using wavelet 
transformation, one can obtain the information of scale and location simultaneously. 

Theory of discrete wavelet transformation is well known as multiresolution approxi- 
mation (MRA) of function spaces. The wavelet bases are remarkably useful because they 
are the orthonormal complete, for some cases unconditional, basis not only of space 
but also of many function spaces such as Sobolev spaces. Holder spaces. Hardy space, 
Besov spaces, etc Dyadic dilation property of the wavelet basis seems quite akin to 
the idea of scaling laws, which appears in many fields of physics such as quantum field 
theory 0, critical phenomena or fully developed turbulence 0. 

The wavelets, however, are scalar function so that application to divergence-free 
vector field contains a problem. Surely one can apply the scalar- valued wavelet transform 
to each component of a divergence-free vector field, say u{x) = ( Ux{x,y, z), Uy{x,y, z), 
Uz{x, y,z) ), and obtain the wavelet spectrum of the field, 

A 

where ipx is scalar wavelet and (*, *) denotes inner product. Each term of the spectrum 

ux{x) = {{ipx, u^)il)x{x), {ipx, Uy)il)x{x), {ipx, Uz)tlJx{x)) (2) 

is not a divergence-free vector field in general. This discrepancy has its root in the 
fact that only two of three components are independent, but on the whole they are 
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dependent each other. (Furthermore they depends also on boundary conditions.) Thus 
the divergence-free vector-valued wavelet function is required for practical purposes. 

Divergence-free vector wavelet bases have been proposed by Battle and Federbush0, 
and Prick and Zimin|^. Battle and Federbush adopted such a constructing way that min- 
imizes the value of integral /(V x A)'^dx under the two constraints, the divergence-free 
condition and an appropriate boundary condition. On the other hand, Prick and Zimin 
proposed such a wavelet that is, roughly speaking, given by the curl of the function given 
by Fourier integral of a step function supported on a spherical shell in the Fourier space. 
The former approach requires the variational calculations when the wavelet transforma- 
tion is carried out. The latter approach, on the other hand, has no such steps because it 
is based on the sharp decomposition of Fourier space into spherical shells. The obtained 
wavelet functions, however, are not orthogonal each other in general, and not localized 
well, in other words, breaks r-regular condition. Because of these discrepancies, the 
approaches they proposed do not seem popular in practical applications. 

In the present work, we propose a general construction procedure of the orthonor- 
mal complete divergence-free vector wavelet basis of L^(T^), in which only two popular 
algorithms, fast Fourier transform (FFT) and fast wavelet transform (FWT), and no 
additional novel one is required. 

The procedure we will describe here is based on a quite different idea from the previous 
two. There are two key ideas for the construction. One is that any function expansion 
using an orthonormal complete basis, {/a; A G A} (say), is a unitary transform from 
L^(T^) to /^(A), where A is an appropriate set of indices. Therefore the Fourier coeffi- 
cients of the base functions fx{k) are regarded as components of an infinite dimensional 
unitary matrix which acts on P space which maps Fourier coefficients to {/A}-expansion 
coefficients. The other one is that all the orthonormal complete basis of the function 
space of solenoidal vector fields on T'^, which is denoted by L|(T'^) hereafter, is given by 
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a certain unitary transform of the complex helical wave basis [0, ^. Thus functions to 
which complex helical waves are unitary transformed by the matrix {J^{fx}} constitutes 
an another orthonormal basis of L|(T^). As an orthonormal complete basis of L^(T^), 
wavelet basis is adopted here. 

This study is an attempt to construct the multiresolution approximation of the vector 
fields. In the present study we restrict our interest on the vector fields on the three-torus 

and the possibility of construction of multiresolution approximation of them. The 
reason of choice of the manifold is mainly due to the fact that the Hodge decompo- 
sition theorem is established on bounded manifolds. Thus we base our attempt on the 
MRA of L^(T'^) though the theory of wavelets is firstly established on the unbounded 
Euclidean space R^. 

This paper is organized as follows. In §2 notations in the paper are explained. The 
orthogonal decomposition of the vector field on is proved in §3. Algorithm of helical 
wavelet decomposition is given in §4, which is an anthology of the properties of helical 
basis. Construction the theory of MRA of vector fields is tried in §5. The Riesz basis 
condition, which is one of the basic properties of MRA, is shown to be broken. Section 
6 is devoted to the discussion on the regularity, ie the localization property of the helical 
wavelet. Finally some remarks are given in §7. 

2 Nomenclature 

Before going into the details, some notations should be fixed. 

Let us denote by R'^ a linear space spanned by a Cartesian basis {e^, Sy, Sz}, and 
the quotient space R^/Z^. x(T^) is a set of vector fields on defined as 

XiT^) ■=lu; u{x)= u^,Uy,Uz e C°°{T^) \, 

i=x,y,z 

where {e^, Cy, e^} is the basis of the tangent space TxT^ at the point x E obtained 
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by the canonical identification. In other words, xC^^) is the set of all the sections 
of the tangent bundle TT^. We identify, if required, the space with the set of 1-forms 
n\T^), or that of 2-forms ^^(t^). 

The completion of F with respect to the norm of the Banach space E is denoted by 
clos£;{F}. In the following analysis, we study the function space given by the L^-norm 
completion of x(T^), 

L|.(T3) := clos {x(T=^)}, (3) 
= |w ; u{x)= Ui{x)ei, u^,Uy,u^ e L'^{T^) >, (4) 

i=x,y,z 

which is a Hilbert space equipped with a inner product, 

{u, v)^ := j^^ u{x) ■ v{x) dx = j^^ (^u.^{x)v^{x) + Uy{x)vy{x) + U;,{x)v^{x)^ dx, (5) 

where u, v E L^(T^), • denotes the scalar product of two vectors, * the complex conju- 
gate, and da; the Lebesgue measure on T^. By definition, any Hilbert space is a Banach 
space with the norm naturally determined by its inner product. Concerning L^(T^), 
the norm is 



||w|L^(T^)|| := sj{u,u)^ = j^^{\u^{x,y,z)\'^ + \uy{x,y,z)\'^ + \u^{x,y,z)\'^)dx, (6) 

where || * |£'|| denotes the norm of a Banach space E. In the following, we drop the 
symbol in definite integrals, and represent by a periodic unit cube, i.e. = [0, 1]^. 
Fourier series representation of f{x) e L^(T^) is formally written as 

/N= E /(fc)exp(27ri/c-a;), (7) 
feezs 

where f{kys {k G Z^) are Fourier coefficients. The calligraphic letter T is used to denote 
the sequence of Fourier coefficients, i.e. = {f{k);k e Z^}. It is also used for the 
Fourier transform of a set of function and a function space, for example ^L^(T^) = 
{.F/;/eL2(T3)}. 



Let us consider a trivial bundle M := x R^. Fourier transform of a vector field is 
defined by the Fourier transform of the components with respect to the Cartesian basis 
as follows: 

u{k) := Ux{k)e.x + Uy{k)ey + Uz{k)ez, (8) 

where k E and {e^, Sy, e^} is a Cartesian basis of the fiber B? of M. Thus (fc, u{k)) e 
M and Tu is a sequence of three dimensional vectors. 
Being introduced the Z^-norm of J^L^(T^) defined by 

V feezs 

J^L^(T^) becomes a Banach space. Applying Parseval identity to each Cartesian com- 
ponent Uy and of a L^(T^) vector field tt, we conclude that the L^(T^) norm 
and ^^(Z^) norm are equivalent. Due to this equivalence of two norms, the fundamental 
sequence of u{x) defined by 

Un{x):^ u(fc)exp(27rife • a;), e N. (10) 

o<|fe|<Ar 

is a Cauchy sequence of vector fields in the sense of L^-norm. In the following, the 
derivatives are formally defined by Fourier series, 

■ \ 

—^■=—— 2mkjUi{k)ex^{2mk- x), (11) 

(JXj l^^J feez3\{o} 

where i, j — x, y and z. 

We distinguish the terms divergence- free and solenoidal in the present work; the 
former is used for such vector fields u that satisfy V • w — 0. The latter term, on 
the other hand, is used when a vector field u is given by curl of certain vector field a: 
w = V X a. Difference of these two kinds of vectors persists in the fact that the harmonic 
vector fields, which are constant function for the case of T^, are also divergence- free. In 
terms of forms, divergence-free fields correspond to closed 2-forms and solenoidal fields 
exact 2-forms. 
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3 orthogonal decomposition of the vector fields on 
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In order to construct complex helical waves, a triad of vector fields {er{k), e^(fc), e^(fc)} 
which is almost identical to a spherical coordinate system basis is introduced in the 
Fourier space. In the present study, they are defined in terms of the wavenumber vector 

k and Cartesian coordinate system basis {6^,6^,62} as follows: 

(erik) X (e,.(fc)/e^) 
e, X eJk) , , ' e^(^) '■= ev(^) X ^r{k). (12) 
--— (otherwise) 
X e.r{k)\ 

The helical vectors hs{k), where the index s denotes polarity of them and is +, — or 0, 
are a triad of complex valued vector fields on the Fourier space that are defined by 

We also use the notations hY:+{k), /is-(fc) and h£){k) instead of h^{k), h_{k) and 
ho{k), respectively. 

It should be noted here that we used definition of helical vectors which is used in 
Ref. with slight modifications, introduction of ho{k) and normalization of h±{k) vec- 
tors. The helical vectors are defined on R'^\{0}, and infinitely differentiable vector fields 
on R^\{0} except on the line along the north pole and the south pole for x and y com- 
ponent of h±{k) and around k = for z components of hj-{k) and all the components 
of hQ{k). In this section, however, we treat the vector field on T'^ and restrict the case 
for k e Z3\{0}. 

Utilizing the helical vectors, we define the complex helical waves that are complex 
valued vector fields by 

V{k,s;x) := hs{k)exp{2TTik ■ x), (14) 

for k G Z^\{0} and s = +, — , 0. According to the polarity of the helical vectors, we 
will call each helical wave S — and D-mode, respectively. It should be remarked 



that we introduced D-mode of complex helical waves that is not defined by Lesieurp] or 
Waleffe[^]. We will see in the following that this mode is curl-free. 
The helical vectors are defined to satisfy the orthonormal relation 

hjk) -hs'ik) =6{s\s'), (15) 

for s, s' = +, — , at each k. Hereafter 6{A\B) denotes Kronecker's delta whose value is 
one only when the two arguments A and B, which are not only numbers but also vectors, 
symbols, etc., coincide, otherwise it gives zero. Therefore the complex helical waves are 



such vector fields that are orthonormal in L^(T^): 

{V{k, s; x), V{k', s'; x))^ = 6{k,s \k' , s') . (16) 

Scalar and vector products of the wavenumber vector and the helical vectors are 

ik-h±{k) = 0, ik ■ ho{k) = \k\. (17) 
ikx h±{k) = ±\k\h±{k), ikxho{k)=0, (18) 

Corresponding to these products, divergence and curl of complex helical waves are 

V ■V{k,±;x) = 0, V ■V{k,0;x) = 27i\k\exp{27iik-x), (19) 

V xV{k,±;x) =±27i\k\V{k,±;x), V x V{k,0; x) = 0. (20) 



The E±-modes are eigenfunctions of the curl operator which are associated with the 
eigenvalues ±27r|fc| and the Z^-mode belongs to the kernel of the curl operator. Thus 
it is obvious now that the indices of helical vectors +, — and correspond to the signs 
of eigenvalues, and that the implication of the alternative indices S+, S— and D are 
solenoidal with positive helicity, solenoidal with negative helicity and dilatational, respec- 
tively. 

Since er{—k) = —er{k), e^( — fc) = e^(fc), e^pi^—k) = — e^(fc), the complex conju- 
gates of helical vectors at k are equal to ones at — fc: 

h^ik) = h,i-k). (21) 
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Therefore, complex conjugates of the complex helical waves satisfy 

V{k,s;x) = V{-k,s;x), (22) 

for k G Z^\{0}, s = +, —, 0. The S+ and S— vectors are complex conjugate each other, 

Mfc) = ^^(fc), (23) 

so that the E±-modes of the complex helical waves are complex conjugate each other: 

Vik,±;x) = Vik,T;x). (24) 

Using these helical vectors, we define the helical decomposition of Fourier coefficients 
(or IP^- sequence), say {u{k)}, by 

u{k) = u+{k)h+{k) + u_{k)h_{k) + uo{k)ho{k), (25) 

where Us{k)^s are the s-mode helical Fourier coefficients defined by the scalar product 

Us{k) ■=u{k) -hjk), (26) 

for each k G Z^\{0} and s = +, — , 0. Each helical Fourier coefficients {us{k)} is a 
sequence of scalars. 

Conversely, by multiplying helical vector on each terms of a scalar (P-)sequence 
{u{k);k G Z'^\{0}}, one can make a vector sequence {u{k)hs{k); k G Z'^\{0}} which 
we call the helical pull up of a sequence {u{k)} to s-mode hereafter. 

It should be remarked here that the helical decomposition of solenoidal vector field 
on has also been discussed by Constantin and Majda[0. In their paper, the de- 
composition is given by u±{k) = u{k) ± ier{k) x u[k) for each k G Z'^\{0} under the 
assumption that the vector field is solenoidal k-u{k) = 0. It is easy to see that, under the 
same assumption, the identity u±{k) = \/2u±{k)h±{k) is satisfied for each k G Z'^\{0}. 
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The helical decomposition of Fourier coefficients is unitary in the sense that the 
identities 

u{k) ■ v{k) = u+{k)v+{k) + u_{k)v^{k) + uo{k)vo{k) (27) 

are satisfied for each k e Z^\{0} if {tt(fc)} and {v{k)} belongs to 1^{Z^\{0}). This 
property leads to the identities 

+ \uy{k)\^ + = \u4k)\^ + + |no(fc)P (28) 

for every k G Z'^\{0}, and to the following consequences. 

Proposition 1 (convergence of /^-sequences) If a three-vector sequence {u{k)} be- 
longs to /^(Z'^), each of its hehcal projections {us{k)} is P {Z^\{0}) -sequence for s = +, 
— and 0. Conversely, if an sequence {u{k)} belongs to /^(Z'^\{0}), each of its helical 
pull up {u{k)hs{k)} is 1'^{Z^\{0}) -sequence for s = +, — and 0. 

Proof: equation (|28| ) leads to the inequalities 

\us{k)\' < \uAk)\' + \uy{k)\' + \u^{k)\\ 

for every k G Z^\{0} and s = +, — and 0. Thus the /^-norm of {us{k)} satisfies 

E < E (l^x.(fc)P + \uy{k)\' + \u.{k)\') < \\u\ll{Z')\\ < oo. 

fc6Z3\{0} fcez3\{0} 

The latter part is proven by two steps. First, the absolute values of the helical vectors 
are one, i.e. \hs{k)\ = 1, so that all the Cartesian components of them are equal to or 
less than one, i.e. \{hs{k))j\ < 1, for every k G Z'^\{0}, s = +, — , 0, and j = x, y, z. 
All the Cartesian component of helical pull up satisfies 

\{u{k)hm,\ = \m\\{Km,\<\n{k)\ 

so that they are /^(Z^\{0})-sequence. Thus the helical pull up is an /^(Z'^yjOD-sequence 
by definition. □ 

Thus the helical transform is unitary transform acts on /^(Z'^\{0}). 
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Theorem 2 (convergence of helical vector fields) If a scalar sequence {u{k)} be- 
longs to /^(Z^\{0}), the vector field that is given by 

Us{x) := ^ u{k)hs{k) ex'p{2nik-x) (29) 
fcez3\{o} 

converges in the sense of L'^(T^)-norni for each polarity s = +, — and 0. 

Proof: Applying Riesz-Fisher theorem to each component, they are proved to belong 
to L^(T^). Thus the vector field belongs to L^(T^) by definition. □ 

We call the vector field expressed by Eq.(p9D s-mode helical vector field. This theorem 
allows us to define the function subspaces of L^(T'^) by 

Ll{T^) ■= J J2 u{k)hs{k) exp (27rifc ■ x) ] {u{k)} G f {z\{0}) i , (30) 
ifcez3\{o} J 

and the projection operators Vs that are linear maps from L^(T'^) to L^(T^) by 

Vsu{x) := {u{k) ■7ijk)^hs{k)exp{2mk- x), (31) 

fceZ3\{o} 

for s = ,D. We call the subspaces helical vector spaces in general, and call each 

one S — , D-mode space, respectively. 

The helical vector spaces L|^(T^), L|_(T^) and L|,(T^) are Hilbert spaces w.r.t. the 
inner product (*)^ and orthogonal each other, that is, the inner product of arbitrary two 
vectors u G L^(T'^) and v G L^,(T'^) is zero if the helical mode indices do not coincide 
s 7^ s'. Therefore the following relation holds; 

^l+(T^) © ^I-(T^) © Ll{T^) C Ll{T^). (32) 

The orthogonal complement of the direct sum in L^(T^) is the space of the harmonic 
functions on T'^. We will discuss it later. This orthogonal relations are expressed in 
terms of the projection operators by 

VsVs' = S{s\s')Vs, (33) 
11 



for s, s' = S+, S— , D. 

The theorem also enable us to define the pull up operator Vl which is a map from 
the space of scalar functions L^(T^) to the space of s-mode helical vector fields L^(T^) 

by 

Vlu{x):= u{k)hsik)exp{2nik-x), (34) 

fcGZ3\{0} 

and their adjoint operators P^* by 

Vl*u{x):= {uik)-hs{k)) exp{2mk-x), (35) 

fcez3\{o} 

which satisfy the identity 

{Vlf, g)^ = (/, rl*9) , (36) 

for every / G L^(T^) and g G L^(T'^). Each pull up operator is unitary in the sense 
that the identity 

{f,9) = m9iO) + {Vlf,Vlg)^ (37) 

holds for every f,g& L^(T'^) and s = +, — , 0. Let us define the zero-mean function 
space Ll{T^) by 

Ll{T^) := i^u{x) ; u{x) G L'^{T^), J u{x)dx = o| . (38) 

The space is Hilbert with respect to the inner product Eq.(^), and the pull up operators 
become the unitary operators from Lq(T^) in the strict sense. This unitary nature of 
the pull up operators allows us to construct an orthonormal complete basis of s-mode 
space from a scalar function basis. 

Theorem 3 (construction of the helical basis) If {fx A G A} is an orthonormal 
complete basis of Lq(T^) where A is appropriate set of indices, then the pull up of the 
basis to s-mode space V^lfx} is an orthonormal complete basis of L^(T^) where s = S+, 
D. 
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Proof: Orthonormality is obvious because of the unitary relation Eg. (pTf) . Complete- 
ness is proved as follows: if {Vlf\,Vlg)^ = for all A G A, then g = 0. Because the 
unitary relation {Vlfx^Vlg)^ = {fx,g) holds and {fx} is complete. Therefore Vlg = 
is concluded. □ 

The theorem leads to an important corollary. 

Corollary 4 (construction of the solenoidal basis) Under the same conditions as 
the previous theorem, the union of the pull up of the basis to S+- and S — mode spaces, 
"^e+I/a} U Px;-{/a}> ■'s cin orthonormal complete basis of the function space of the square 
integrable solenoidal vector fields. 

Finally, we shall discuss the "residual" of the projection operators in order to salvage 
the modes with wavenumber 0, which are left behind the helical modes. It is easy to see 
that the function defined by 

Uhix) := u{x) - + Ps- + Vd)u{x) (39) 

is a uniform vector field on for every u{x) G L^(T^). Because 

u^{x) - + Pe_ + Vd)un{x) = M,(0)e, + Uy{0)ey + M,(0)e, (40) 

for every elements of the fundamental sequence of u{x). We will denote the space of 
uniform vector fields by and define the projection operator Vh from to by 

Vhu{x) := u^{0)e^ + Uy{0)ey + ^2(0)62. (41) 

By definition, ^^(T^) n ^^(T^) = {0}, for s = S+, S-, D. Every elements of ^^^(T^) 
is a harmonic function on T^, i.e. satisfies Vh = 0,Vxh = O. The degree of freedom 
of L'jj{T^) is three. It is known that the number coincides with the first Betti number of 
the manifold T^p|. 

At this point, the orthogonal decomposition of L^(T^) is completed. 
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Theorem 5 (Hodge-Beltrami decomposition) 

L^(T^) = L|+(T^) e Ll_(T^) e Ll(T^) Ljj(T^). (42) 

This is the synthesis of a special case of Hodge decomposition and the Beltrami decom- 
position. 

Let us finish this section by giving a theorem to construct the orthonormal complete 
basis of divergence-free vector fields. For this purpose, we introduce a triad of uniform 
vector fields { /is+(0), /is_(0), /id(0) }. Choice of the vector fields is arbitrary except 
for the requirement that orthonormal relations hs{0) ■ hs'{0) = S{s\s') must be hold. 

Theorem 6 (construction of the divergence-free basis) If {fx} is an orthonormal 
complete basis of L^(T^), then 

{h(0)h^+(0)+Vl+fx{x)} U {hiO)h^40)+Vl_fx(x)} U {hnm (43) 

is an orthonormal complete basis of divergence- free vector fields L|^(T^) © L|_(T^) © 

4 helical basis 

Here we gives an anthology of the properties which are satisfied for every hehcal basis. 

First, we describe the procedure to obtain the expansion coefficient of the helical basis. 
Consider the helical basis which is obtained by the helical pull up of an orthonormal basis 
of Lo(T^), say wavelet basis {^a; A G A} where A is an appropriate set of indices. The 
expansion coefficients of an L^-vector field u{x) with respect to Vl{ipx}, the helical pull 
up of {■0a} ^-re given by the inner products 

ux,s := {Vl^Jxix), u{x))^ 
= (V'A(a;),nM^)) 

= |A(fc)(w(fc)-Mfe)), (44) 

fcez3\{o} 
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for A G A. The last equation leads to the following procedures; (1) calculate jFw, the 
vector- valued Fourier transform of u{x): (2) calculate the inner product u hg, and we 
obtain s-mode hehcal Fourier coefficients {us}: (3) calculate J^~^{us}, the scalar-valued 
inverse Fourier transform of the s-mode helical Fourier coefficients: (4) applying the pro- 
cedure of the scalar-valued function expansion, for example, the fast wavelet transform 
algorithm to J^''^{us}. Thus we obtain the expansion coefficients {ux,s}- It should be 
remarked here that the procedure given above requires no novel and specific algorithm, 
but only such tools as FFT and FWT, which are popular in the signal processing. 

The helical pull up of a real valued scalar function f{x) e L^(T^) is also a real valued 
vector field. Because the helical vector satisfies 

h,{k) = hs{-k) (45) 

for each fc e so that each component of the helical pull up of the Fourier coefficient 
f{k) satisfies 

(mhsik)), = [f(-k)hs(-k)) , (46) 

for j = 1, 2, 3 and each k e Z^\{0}. Thus each element of the hehcal basis pulled up 
from a real- valued function basis of L^(T^) is also real- valued. 

The continuity and differentiability of each element of a helical basis depends on those 
of the original scalar function basis. These properties of the helical basis, however, may 
not be identical to those of original one. Since the multiplication of the polarity vector 
hs{k) lead to the 'twisting' of the values of the scalar function to each components. 
In order to evaluate the properties, here we go via such a path that goes through the 
Sobolev's imbedding theorem. 

When the function f(x) belongs to the Sobolev space H^{T^), i.e. f{x) satisfies 

\\f{x)\ff^-{T')\\'= Y.i^ + m^m\'<oc, (47) 
feez3 
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each component of its helical pull up satisfies 

\\{Vlf{x))^\H%T')\f= E (1 + \k\rmk)\'\ihs{k)h\' < \\f{x)\H^{T'W (48) 

feez3 

for j = X, y and z. This convergence property and the Sobolev's embedding theorem 
guarantees the continuity of the helical pull up of f{x) when s > |. For this case, the 
formula of the helical pull up Vlf{x) converges pointwise. The differentiabihty of the 
helical pull up of f{x) is guaranteed up to r-th order when f[x) G H^^i(T^) where 
s >r. 

Next we consider the relation between the E+-mode and the E — mode associated 
with exchange of parity, ie reversion of orientation of the coordinate system on T"^. The 
parity exchange operator is denoted by F here. When F acts on a real-valued scalar 
function, say u{x), the relation 

ru{x) = ^ u{k) exp(27rife • {-x)) = ^ u{k) exp(27rife • x) (49) 
feez3\{o} feez3\{o} 

holds. Thus the mirror image of a real-valued scalar function u is given by the inverse 

Fourier transformation of the complex conjugate of the Fourier transform of u, that is 

Tu{x) =J^-^ (Tu){x). (50) 

The mirror image of the helical pull up of a real- valued scalar function is given by 

TVl^u{x) = M(fc)/iE±(fe)exp(27rifc- (-x)), (51) 

fe6Z3\{0} 

= ^ u{k)h^^{k) exp(27rifc • x). (52) 
feez3\{o} 

This leads to the commutation relation between the parity exchange and the helical pull 
up for real- valued functions: 

TT'L = ^Lr. (53) 
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5 mult iresolut ion approximation of the vector fields 

In this section, the properties of the MRA of Lq(T^), which is defined by the restriction 
of the MRA of L^(T^) to Lq(T^), and its hehcal pull up's are discussed. The omission 
of constant functions from the MRA causes the absence of such solenoidal basis that is 
homogeneous, i.e. given by the orbit of a function by the action of finite group. It is 
discussed that this "discrepancy" is consistent with the practically natural postulation 
for approximating the constant vector field on the finite number of grid points. 

First, we briefly review the multiresolution approximation (MRA) of the spaces of 
the functions of period 1. For details, one should consult §3.11 of Ref.[|I[]. The MRA of 
the function spaces on T is obtained by the periodification of the wavelet functions. 

Consider the r-regular MRA of L^(R), denoted by {V^(R); j G Z}. The completion 
of each subspaces V^(R) with respect to the L°°-norm, cIosloo(j^) {Vj (R)}, retains the 
relation /(x) G Vj <^=^ /(2x) G Vj+i, which is one of the MRA conditions. Let us define 
V,{T) by 

VjiV ■■= fix) G ^clos^ {VAm, fix + 1) = . (54) 

The nested sequence {Vj{T)] j G Z} is called the r-regular multiresolution approximation 
ofL\T). 

According to the lemma 13 in §3.11 of Ref.[]T|, the spaces V^(T) have the following 
properties. If j < 0, they are identical. The space Vo(T) consists of constant functions. 
The dimension of the space ^ (T) is 2-' . 

Let Wj{T) be the orthogonal complement of Vj{T) in V^+i(T), the space L^(T) is 
represented in terms of a direct sum of the subspaces as follows: 

L\T) = VoiT) © WoiT) © l^i(T) © W2iT) © .... (55) 

In the present study, we will define the multiresolution approximation of L^(T^) by the 
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tensor product of the one-dimensional ones: 

Vj{T^) ■.= Vj{T)^Vj{T)^Vj{T) for J = 0,1,2,.... (56) 

The corresponding spaces Wj{T^) {j > 0) arc defined by the orthogonal complement of 
V^(T^) in V^_|_i(T^). According to the construction procedures, Vb(T^) consists of the 
constant functions. The dimension of T^ (T^) and Wj{T^) are 2^-' and 7-2^\ respectively. 

Now we will try to swim away from the shore of the established MRA theory in order 
to find out the hidden bank of the MRA theory of the three-dimensional vector field. 

In section three, we introduced the zero-mean function space Lq(T^) in order to 
construct the solenoidal function space by helical pull up. In terms of the MRA of L^(T'^), 
one can obtain a nested sequence of the zero-mean subspaces of L^(T^), {Vj{T^)}, defined 

by 

^i(T') = 0M/.(T-^) forj = 1,2,3,..., (57) 

i=Q 

or equivalently defined by the orthogonal complement o/Vo(T^) in Vj{T^). The obtained 
sequence retains the following conditions of MRA by definition: 



^.■(T^') C V;-,i(T=') for^j>l; (58) 
fix) e Vj{T') ^ f{2x) e V;+i(T=^) for > 1; (59) 




It is easy to see that the dimension of l^(T^) is 2^^ — 1. The last condition is easy 
to prove. Because all the functions in the space are, by definition, orthogonal to all 
the constant functions that belong to Vo{T^), that is, Jrj,3 f{x)dx — for ^/(x) e 
cloSi2(T3) {0°li t/,(T^)|. Hereafter we will call the nested sequence {^^(T^)} the MRA 
ofLl{T^) or the zero-mean MRA of L'^{T^). 

In the previous section, we proved that the helical pull up operators are unitary in 
the strict sense when they act on Lq(T^), instead of L^(T^). Thus the nested sequences 
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that are given by the hehcal pull up of the zero-mean MRA of L^(T^), {7^]V^(T'^)}, 
satisfies the conditions 

VlVj{T')cVlV,+,{T')cLl{T') for^j>l; (61) 
u{x) e VlVj{T') ^ u{2x) e VlVj^T') for > 1; (62) 

4(T3) {®^=' ^^^^^^^'^} = (63) 
for s = S-I-, E— and D. We will call each nested sequence the MRA of Ll{T^) or the 
s-mode helical MRA. Because the pull up operator is unitary, the dimension of P]V^(T^) 
is 2^^' - 1 for s = E+, E- and D. 

Intentionally we do not discuss the Riesz basis condition, which is one of the properties 
that constitutes the definition of multiresolution approximation, and r-regular condition 
till now. On the latter condition we will discuss in the next section. 

In order to discuss the former condition, we define the term homogeneous basis. 
Consider a Hilbert space on T"^, say H, and the finite group of residues k/N modulo 
1, Fjv, where m, N & N and k e Z"*. We will say H has homogeneous basis if there 
is such a natural number N and a function 0(x) whose orbit under the action of Fpf, 
{4>{x — k/N) ; k e Z™} is orthonormal complete basis of H. For example, each ^(T) 
has homogeneous basis which is given by the orbit of the scaling function 2^^'^(f){2^x) by 

The dimension of each subspace of zero-mean MRA of L^(T^) Vj(T^) is 2^^ — 1. 
Because there is no such an integer p that satisfies 2^^ — 1 — p^, Vj{T^) never has any 
homogeneous basis. So is its helical pull up VlVj{T^) 

Proposition 7 (absence of the homogeneous beisis) Each subspace of the MRA of 
L^(T^), VlVj{T''^) (s = E+, S— , D and j eN), cannot have any homogeneous basis. 

Intuitively this is a consequence of a practical requirement for the signal processing 
or the numerical analysis that the vector field u{x) should be regarded as an element 
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of L^(T^) if all of its values on the grid points located at x = are identical. In 
mathematical words, if there exists a vector- valued homogeneous basis {(f)ji{x)}, where 
j and I are indices for the resolution and the location, respectively, it should be able 
to approximate not only the elements of L^(T^) but also those of L^(T^). Thus the 
requirement is restated as the following postulation. 

Postulation 8 (homogeneous approximation) If u{x) is approximated as Uj{x) — 
^ Cj4>ji{x), where cj is a constant for each finite resolution class j, i.e. the inner products 
Cj^i = {4>jii u)^ are dependent only on j and independent of I, u belongs to L^(T^). 

Defining the hehcal pull up of the scahng function of j-th resolution (f)j{x) by 

cl>j^,{x):^mhs{0)+Vl(l>{x), (64) 

where ^1.5(0) is an arbitrary uniform vector field with amplitude one, one can obtain the 
set of function that is given by {r2j 0j ^(cc)}. It approximates any constant vector field 
in such a way that is stated above. 

6 singularity of Edz-mode helical vectors and spatial 
coherence of solenoidal helical wavelet 

In this section the spatial coherence of hehcal wavelet is discussed. 

MRA is called r-regular when the generating scaling function ■0o and the associated 
mother wavelets ■^g satisfy the conditions 

Vm e Z, 3Cm < 00 s.t. |9"'0e(^)| < Cm{l + \x\y'^ (65) 

for every multi- index a satisfying |ci;| < r. Our interest here is whether the regularity of 
the helical wavelet suffers from the singularities of the E±-mode helical vectors. For this 
purpose, in this section we will discuss spatial coherence in terms of the vector valued 
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function on defined by the Fourier integral of L^(R^) function: 

^ (x) = / (J^iljx)(k)hs(k)exp(2mk-x)dk. (66) 

We will call the integral transform integral helical pull up (to s-mode) hereafter. Vector 
fields on T'^, which we have discussed, are given by periodification of the function. We 
will assume {J^il'x){^) = in order to avoid the arbitrariness of the definition at fc = 0. 

Our afraid is as follows. In order to construct solenoidal basis by hehcal pull up from 
an orthonormal complete scalar function basis {ipx] A G A}, the conditions 

k ■ ^s±(fc) = 0, (67) 
|^s±(fc)| = l, (68) 
J2\J^Mk)f = l, (69) 

A6A 

are required for every k. The first condition says that the vector field is tangential to 
spheres = const.. It is a well known result of the differential topology that two-sphere 
S"^ is not parallelizable manifold so that the vector fields on S"^ have at least two singular 
points]!^. Thus h-£±{k) must have singular points for every |fc|'s. The second condition 
says that the vector fields should behave like the function sgn(a;) around the singular 
points. The third condition together with the two previous ones leads to a conclusion that 
there must exist such functions that have singular points in their supports. Therefore, 
even if ipx belongs to the Schwartz class iS(R^), the vector function J^4'xh-£± breaks 
the continuity and does not belong to S when it has singular points in its support. 
Furthermore, these properties leads to an expectation that the inverse Fourier transform 
of J-'ipxhY.± has algebraic tail like the Fourier transform of sgn function. 

Asymptotic analysis of the Fourier transform of S±-mode helical vector at large 
r shows the far field behaviour as around 2; ~ plane (see Appendix A). The 
same exponent of algebraic tail is obtained by a simple scaling argument and the two- 
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dimensional distribution of the tail is required by the convergence of wavelet function 
under the periodification operation (see Appendix B). 

Thus it is plausible that, in the orthonormal solenoidal basis on which is con- 
structed by the integral helical pull up of a scalar basis, there exist such functions that 
have algebraic tail. Thus we conjecture that there exist no such orthonormal solenoidal 
wavelet basis that all the species of wavelets decreases rapidly at large r. 

In the preceding part of this section we will examine the behaviour at large r numer- 
ically using Littlewood-Paley type MRA and its helical pull up. 



According to the recipe by Yamada and Ohkitani [13|, the Fourier image (f){k) of the 
scaling function (j){x) is given by 



(70) 



where 



h{\-k) 



h{k) 



exp(-A;"2) /c > 0, 
k<Q. 



(71) 



h{k ~ \) + h{i - ky 

By definition, this scaling function is infinitely differentiable and following properties are 
satisfied: 



for - i < A; < i, 
< 0(/c) < 1 for - I < < -i, \<k <\ 



m = 1 

< ^ 
${k) = 



3 " " 3' 3 

for < -|, l<k, 



2 

3' 



and 



5:i0(fc+j) =1. 

Then, the mother wavelet ipi^) is obtained by the Fourier transform of the image, 



(72) 



(73) 



\ 



'k' 



{k f exp(— Tri/c). 



(74) 



Three dimensional wavelet functions are constructed by tensor products of one dimen- 
sional ones. There are seven species of mother wavelets i/'e (e = 1, 2, 7) given by 



xlj,{x,y,z) = ip^{x)ipn{y)i^ciz), (^, i], C = or 1), 



(75) 
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where ijJo{x) '■= (pix) is the one dimensional scahng function, ipi{x) := ipi^x) the one 
dimensional mother wavelet, and the label e is determined as e = ^ + 2r^ + 4C. 

It is obvious from the definition given above that Meyer's wavelet belongs to Schwarz 
class. Because its Fourier transform has compact support, are bounded and infinitely 
differentiable. Thus the Littlewood-Paley type MRA (Meyer's wavelet) is r-regular for 
any positive integer r. 

According to our definition of helical vectors, x and y components of type e = 4 
helical Meyer wavelets Vl2±ipo{x)ipo{y)'ipi{z) have singular points in their support and do 
not belong to Schwarz class. Spatial coherence of these two components are not obvious. 
It should be remarked that the Fourier transform of these components remain bounded 
and compact supported so that Sobolev's imbedding theorem guarantees the analyticity 
of them. All the components of helical mother wavelets except these two retain the 
r-regular conditions. 

In order to evaluate the localization of wavelet functions, we introduce coherence 
spectrum defined by 



where the origin of spherical coordinate system (r, (p) is taken at c. 

In Fig.l, the coherence spectra of the type-e wavelet functions are depicted. The other 
parameters that define the wavelet are the resolution class j, location I, and helicity s. 
Those of the tested wavelet are set to j = 6, Z = and s — +, respectively. The center 
of the spectrum c is taken at (2"''"-'^^, 2~^~^?7,2"''~^C) for the type-e scalar or hehcal 
mother wavelet, where scalar Meyer wavelet ipeeoi^) takes its maximum. The function 
is evaluated numerically on the 256^ number of grids given by a; = {ix/N,iy/N,iz/N) 
where N — 256, ix, iy and are integers that satisfy —N/2 < ix-,iy,iz < N/2. The 




(76) 
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spectrum is approximated by the sum 




(77) 



rj-A<\X-C\<rj+A 



where the radius Vj and the shell thickness A are rj = j/N and A = 1/2N, respectively. 

In Fig. 1(a) the spectra of scalar wavelet functions are depicted. All the wavelets 
show rather exponential behaviour at large r with oscillations. Helical wavelets, except 
for the type-4 one, show almost the same behaviour as the scalar ones (see Fig. 1(b)). 
The oscillations at large r are shghtly modified. There are two remarkable features in 
the spectrum of the type-4 helical wavelet. One is that the spectrum decreases more 
slowly than those of the other wavelets. The functional form of the spectrum seems 
rather algebraic at medium scales about 0.05 < r < 0.3. The other is that the relative 
amphtude of the oscillations is smaller than that of other hehcal wavelets. In Fig. 1(c), 
the contribution of each component of type-4 wavelet to the spectrum is depicted. It 
is obvious from the figure that x- and y-components contribute to the features given 
above. The 2;-component, on the other hand, shows the same rapidly decreasing features 
as others. 

On the whole, the type-4 hehcal wavelet is less locahzed and less oscillating at large r 
than the others. These numerically investigated features are consistent with the results 
of asymptotic analysis and scaling argument given in the appendices. 

7 concluding remarks 

In the present paper, we proved the Hodge- Beltrami decomposition of vector fields on T^. 
Then, we have shown the construction procedure of the orthonormal complete basis of 
the solenoidal fields on T^. This procedure has a merit that it requires only conventional 
numerical algorithms (ie fast Fourier transform and fast wavelet transform) and no novel 
one in calculating the wavelet coefficients. 
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Based on this procedure, multiresolution approximation (MRA) of vector fields on 
is constructed. It was shown that there exists no scahng function which has close 
relation to the assumption that if a vector field is constant in numerical simulation, it 
can be regarded as an approximation of a uniform field. Nonexistence of a single scaling 
function which generates MRA might have relation to the idea of multiwavelet . 

We have conjectured that, in any divergence-free three-dimensional vector wavelet 
basis, there must be such kind of wavelet function that behaves algebraically at large 
I a; I so that it breaks the r- regular condition. This conjecture has its root in the facts 
that the helical vector must be tangential to 5*^ for every point, and that 5*^ is not a 
parallelizable manifold. We feel that the requirement of exponential decay imposed in 
the work by Battle and Federbush is inappropriate if it is required for all the functions 
that constitutes a solenoidal basis . In two or four dimensional spaces, existence of the 
r-regular solenoidal vector wavelet basis is plausible because and are parallelizable. 
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Appendix A: asymptotic analysis of Fourier transform 
of helical vectors 



In this appendix Fourier transform of the hehcal vectors is considered. Since the hehcal 
vectors have singularities, there exists such type of hehcal wavelet that does not belong 
to Schwarz class. In order to estimate the amplitude of helical wavelet at large r, first 
we calculate formally the Fourier transform of each component, and then, asymptotic 
analysis at large r is carried out. 

In Cartesian coordinates, the basis vectors of the spherical polar coordinate system 
are decomposed as 

, , p cos 13 p sin j3 ( 

eJk) - 4^2^ e + 4i£=e - ^ r (78) 

e-^{k) = —smPex+ cos(3ey. 

where (p, /3, Q denotes cylindrical coordinates. We consider here the x-component of the 
Fourier transform of hj]j-{k). By operating the ninety degree rotation around 2;-axis, 
one can obtain the y-component of the Fourier transform. These two components are 
relevant to the singularity problem of wavelet function. 

The Fourier transform of the x-component of e^(fc) is given by 

(e^(fc)) exp (27rifc • x) dk 

/oo f-2iT POO COS /3 / \ 

J (^2i J —^^===^dCj[2cosa j cos/?exp(27ripr cos/?)d/?jpdp 

-AiTcosa p'^Ki{2TTzp) Ji{2'iTpr)dp, (79) 
Jo 



where we define the value at —z by exchanging the sign of Eg . ([79D because the integrand 
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is an odd function. Similarly, the Fourier image of the x-component of e^p^k) is 



R3 



[e^{k))^ exp (27rifc ■ x) dk 







/ / sin/5exp(27ri(pr cos(/3 — a) + C-z) jpdp d/3d^ 

exp(27riCz)dCj f y sin(/3 + a) exp(27ripr cos /3)d/5 ) p dp 



/■oo 

= -27ri5(z) sina / p Ji(27rpr)dp. (80) 

JO 

/^From Eq.(|79D and Eg. (pOf) , one formally obtains the Fourier transform of {h±{k,))x as 
/ ih^ik)) exp (27rifc ■ x) dfc 

= V27r /°?-2piri(27r;zp) cosa±5(;z) sina)pJi(27rpr)dp. (81) 

JO 

Because Dirac's delta function 5(z) is zero for \z\ > and asymptotic behaviour of 
modified Bessel's function Ki{2iTpz) ~ exp(— 27rpz) for \z\ ^ 1, function is well localized 
in the 2;-direction. So we focus on the behaviour around \z\ <^ 1. In the region 271 pz <^ 1 
and 27rpr ^ 1, the Bessel functions are approximated as 

1 



Ki{2Tipz) 



2Tipz 



, Ji(27rpr) 



1 Stt 

cos 27rpr 

vr^/pr V 4 



(82) 



Substituting this expression, we obtain the asymptotic expression of Eq.(|8TD given by 

J (/i±(fc))^ exp (27rifc ■ cc) dfc 

~ ]J- ^ p.v. ^-^ cos a ± S{z) sina^ J cos ^27rpr — dp. (83) 

Let us investigate the behaviour of the definite integral in Eq.(p3D by replacing the semi- 
infinite interval [0, oo) by the finite one [0, po], i.e. 

rpo f 37r\ 

/(po) := cos \2'Kpr - — j dp. (84) 

It is easily confirmed that 

/(po) = J^Vp cos i^27Tpr - Y j dp 



2nr 



Po sin (27rpor - 



4 / 2^2^ 



S{2^) + C{2^) 



^5) 
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where S and C are Fresnel integrals, 



cos i^t^ ] dt. 
V2 



For sufficiently large x, they are asymptotically evaluated as follows: 



S{x),C{x)r^^ + o(^ 



Substituting Eq.(^7|) into Eq.(^5l), one obtains 



— v^sm [27ipor- — 



1 



1 1 

T + — 



which leads to the following approximation: 



Po sin 27rpor 



2V2r 



1 + 



1 



^6) 



(87) 



R3 



[h,±(fc))^exp (27rifc ■ cc) dfc 



- { p.v. {- \ cosa±5(z) sin a 

r \ Tx \z 

1 



X lim 



— " T H 

po^oo [4y/27rr2 27rr 



Vp^sin ( 27rpor - + ^ ("^ 



The term 



4'7rr^ V TT 



-p.v. {- ) cos a ± sin a 



(89) 



(90) 



in Eq.(|89D does not depend on po- Thus we can draw a conclusion that the Fourier 
transform of {h±{k))x is a function asymptotically behaves like (x^ + y^)^^ when z ~ 
and I a? I ^ 1. 



30 



Appendix B: some estimations a priori 

In Appendix A we showed the corroboration of the existence of algebraic tail whose 
exponent is —2. The tail is shown to distribute rather two-dimensionally around z ~ 
0. In this appendix we will show that the same conclusion is obtained by a simple 
scaling argument and by the requirement of convergence of wavelet function under the 
periodification operation. 

Assuming that singularities of the helical vectors cause algebraic tail of helical wavelet 
at sufficiently large r = \x\, exponent of the tail is determined by a simple scaling 
argument. When wavelet function in Wj is assumed to be divided into two parts, rapidly 
decreasing part and algebraic tail, as 

Mx)-fix)+Ar-^, (91) 
for sufficiently large r, the wavelet in Wj+i must be given by 

V¥'ijjx(2x) = V¥f(2x) + 25-" Ar-". (92) 

The domain of Fourier integral is enlarged by two in each direction so that the contri- 
bution of the singularities of hehcal vectors to class Wj+i wavelet become twice of those 
in class Wj. (Remember that functional form of h± around the singular points does 
not depend on |fc| and the points are aligned with a one-dimensional manifold). The 
amplitude of the Fourier transform of ijj{2x) is one-eighth of that of tpix). Thus the 
factor of algebraic term should be 2~^, one-fourth of the factor of /, and the exponent 
q; = 2. 

Scaling argument leads to a conclusion that the exponent of algebraic tail is —2. Is 

2 

such behaviour of r isotropically distributed, i.e. found in all the directions of three- 
dimensional space? The answer is no. They distributes two-dimensionally because of 
the following two reasons. One is that the singularity distributes one-dimensionally so 
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that the Fourier transform in the directions perpendicular to the singularity suffers from 
it. The other is that helical wavelet on R^, say tl)ji{x), should not diverge under the 
periodification operation to obtain the wavelet on T^. This requirement leads to an 
estimation of the contribution of algebraic tails from far field to the periodification sum: 

— r^-^ dr < oo ^ D < 2, (93) 

w>ro ° 

where Tq is a sufficiently large number and D is the (fractal) dimension of region where 
the algebraic tails spread over. Though the estimation is quite rough, this result claims 
that the tail must not spread over to all the directions of three-dimensional space, but 
should be confined in the region which is at most two-dimensional. 
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Figure captions 

Fig.l Coherence spectra of wavelet functions; (a) scalar Meyer wavelet, (b) helical Meyer 
wavelet. Type of wavelet is distinguished by the line; solid line: e — 1, 2, broken 
hue; e = 3, 5 and 6, dashed hne: e = 4, dash-dotted line: e — 7. (c) Spectra of 
each component of type-4 helical wavelet. Solid line: x-componcnt, dash-dotted line: y- 
component, broken line: 2;-component. Solid and dash-dotted lines are indistinguishable. 
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Figure 1(a) 
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Figure 1(b) 
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Figure 1(c) 
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